Unusual Length Dependence of Conductance of Some Oligomers 
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Recent experiment found a quantum length dependence of oligothiophene molecule conductance 
at low bias [Xu et al., nano Lett. 5, 1491 (2005)], the long molecule has large conductance. By 
means of a first-principles method we obtain both the quantum length dependence of conductance at 
low bias and the classical length dependence of conductance at high bias region for oligothiophene. 
In between there is an oscillated conductance behavior. The transport behaviors are determined 
by the distinct electronic structures of the molecular compounds. The various conductance length 
dependence may appear for the organic compounds. Our further investigation finds that the classical 
conductance length dependence in polyphenanthrene dithiolates and another unusual conductance 
length dependence in polyacene ditholates: the quantum length dependence of conductance is at 
the high bias and the classical length dependence of conductance is at the low bias. 

PACS numbers: 73.23.-b, 85.65.+h, 31.15.Ar 



Recently, much interest has been focused on molec- 
ular electronics which is in expectation of new genera- 
tion of electronic devices as so-called molecular devices. 
The fabrications of molecular devices require to under- 
stand transport properties of single molecules, and there 
are much experimental and theoretical works reported in 
these years [HI2I1I1I1IS0,I1I3 Some interesting 
effects attracted much attention: gate effect [lfj, [ill LL2J j 
negative differential resistance 0, 0], and diode rec- 
tifier etc ^3- These effects are expected as the basis of 
molecule device design and the investigation of the molec- 
ular devices are in a rapid progress. Except for these ef- 
fects, the length dependence of conductance of molecule 
wires is also discussed widely [3 E E3 El E E3 • 

Lang and Avouris |l7^ studied carbon-atom wires with 
the first-principles method and found a conductance os- 
cillation with respect to the carbon atom number. Simi- 
lar full ab initio computation of carbon-atom wires per- 
formed by Larade et al. pronounced a result con- 
sistent with that of Lang and Avouris. Analogous phe- 
nomenon was also observed in the metallic Na atom wires 
|23| . Except for atom wires, the oligomer is another in- 
teresting molecular wire. For these oligomer wires, early 
theory presented a semi-empirical law called the conduc- 
tance exponential law, G — G m e~ lL , here L is the length 
of the oligomer, 7 is a damping factor [J^. This law is 
valid for the low bias voltage. Latterly, a full ab initio 
computation was applied to some kinds of oligomers and 
reached results conformed with this law |l9|, |2fJ . By us- 
ing the Landauer formulation combined with the density 
functional theory (DFT) Tada et al. [2l| studied sev- 
eral oligomer compounds and declared that the conduc- 
tances for most of them conform with the exponential 
law, but there are two types of oligomers that disobey 
the law, polyphenanthrene dithiolates(PPh(n)DTs) and 
polyacene dithiolates(PA(n)DTs), with a conductance 
oscillation, similar to that in atom wires. The recent 



experiment of the oligothiophene dithiolates((n)TlDTs) 
shows an unusual length dependence of conductance for 
the bias lower than 1.2V. Although 4T1DT is longer 
than 3T1DT, the conductance of 4T1DT is higher than 
3TlDT's The motivation of our study is to discover 
and understand the unusual length dependence of con- 
ductance for oligothiophene and other organic molecules. 

In this letter we employ a full ab initio method to cal- 
culate the transport behavior for three kinds of oligomers, 
(n)TlDTs, PPh(n)DTs and PA(n)DTs (see Fig.l). Our 
numerical results confirm the non-classical conductance 
length dependence of (n)TlDTs at low bias, which was 
shown in the experiment |l2j . and find a classical con- 
ductance length dependence coexisted at high bias. And 
the results of PPh(n)DTs and PA(n)DTs are far differ- 
ent from the semi-empirical calculation |2lj |. Under low 
bias, both PPh(n)DTs and PA(n)DTs obey the expo- 
nential law qualitatively. With bias expanded into the 
high regime, PPh(n)DTs' length dependence retains the 
classical behavior, the longer the molecule is, the lower 
the conductance is. But for PA(n)DTs, the exponential 
law is broken down under high bias. In the high bias 
range the conductance will be inversely proportional to 
the molecule length, the long molecule wires have higher 
conductance than the short ones, which is opposite to 
the classical length dependence of conducting wires. Our 
computations show that coupling between the molecule 
and Au electrodes plays an important role. The coupling 
could shift and expand the molecule orbitals dramatically 
and change the HOMO-LUMO gap by different manners 
due to their distinct electronic structures. Coupled with 
Au leads, the HOMO-LUMO gap of PPh(n)DTs almost 
keeps a fixed value with the molecule length increased. 
On the contrary, the HOMO-LUMO gap of (n)TlDTs 
and PA(n)DTs will decrease obviously with increase in 
the molecule length. The distance between the Fermi 
level and HOMO or LUMO levels determines the conduc- 
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tance gap of the molecules and thus affects their length 
dependence of conductance. This point will be discussed 
latter in more details. 

After considerable effort was made by many people, 
theories based on the first-principles calculation have 
been built commendably. Here, we combine the den- 
sity functional theory and nonequilibrium Green's func- 
tion theory by means of the self-consistent procedure to 
solve the transport problem of the electrode-molecule- 
electrode open system 0, H, El H El • The advantage of 
such method is its theoretical rigor assured by the self- 
consistent procedure and its sensitivity to bonding and 
surface chemistry at contacts. The NEGF equations are 
as follows 

G R = (E+S-F-Z R (1) 
Si = W-/, (i = l,2) (2) 

p = J dE[G R (f 1 T 1 + / 2 r 2 )G A /2^], (3) 

where S and F are the overlap matrix and Fock matrix 
of the molecular part. Y. R are the left (right) self 

energy. In our computation, the self-consistent procedure 
in Gaussian03 [24| is extended from an isolated molecule 
to the lead-molecule-lead open system. The calculations 
for all the quantities in the open system, e.g., the surface 
Green's functions gi, the self-energies £j, the broadening 
function Ti, and the density matrix p are included in the 
main program of Gaussian03 as the subroutine programs. 
After this procedure is converged, transmission T and 
current / can be obtained 

T = Trace(T 1 GT 2 G^), (4) 

/OO 
dET(E,V)(h(E)-f 2 (E)). (5) 
-oc 

The computation detail is described in where the 
DFT method with B3PW91 and LANL2DZ basis are 
adopted. The Fermi level of Au leads is set as —5.1 eV. 

The calculated I-V characteristics of (n)TlDTs, 
PPh(n)DTs and PA(n)DTs are shown in Fig.l, n is the 
number of thiophene ring (a) or benzene ring (b, c). 
(n)TlDTs with small conductance gap have bigger cur- 
rent than other two oligomers. Its length dependence 
of conductance is unusually non-classical under low bias, 
which is consistent well with the experiment reported 
by Xu et al. recently They measured 3T1DT and 
4T1DT in the bias regime from V to 1 V and found 
that the conductance of 3T1DT is smaller than 4TlDTs'. 
Not only does our calculation agree with their experi- 
ment result, but also it shows a very interesting return 
to the classical length dependence of conductance at high 
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FIG. 1: I-V characteristics of the three oligomers. (a) 
(n)TlDTs show the unusual length dependence in the low 
bias regime and the classical length dependence at bias V = 3 
V. (b) PPh(n)DTs show a classical length dependence. The 
inset illustrates that they obey the Magoga's exponential law 
at low bias V = 0.4 V. (c) PPh(n)DTs present the classical 
length dependence at the low bias regime and the quantum 
length dependence at bias V = 2 V. The inset illustrates that 
they obey the Magoga's exponential law at low bias V = 0.4 
V. 



bias, near V = 3 V. PPh(n)DTs illustrate the classical 
length dependence under the whole voltage range from 
V = 1 V to 2 V [Fig. 1(b)], and the exponential law 
(G = G m e~ lL ) qualitatively under low bias (inset of 
Fig. 1(b)). Fig. 1(c) shows a different conductance length 
dependence for PA(n)DTs. Under small bias voltage, 
the conductance exponential law is still satisfied (see in- 
set of Fig. 1(c)), but at high bias, the exponential law 
breaks down, some longer molecules have higher conduc- 
tance than those shorter ones, which leads to an irregular 
conductance oscillation with respect to the n number. If 
bias increases continually to approach a special value, the 
classical length dependence of conductance will be com- 
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FIG. 2: Transmission spectrum of the three oligomers under 
bias 0.2V (shifted vertically for visibility), (a) (n)TlDTs, the 
HOMO-LUMO gap, which is small, and the distance between 
HOMO and Fermi energy decrease with the rise of n number, 
(b) PPh(n)DTs, the HOMO-LUMO gap, which is large, and 
the distance between LUMO and Fermi energy increase with 
the rise of n number, (c) PA(n)DTs, the HOMO-LUMO gap 
and the distance between LUMO and Fermi energy decrease 
with the rise of n number. 



pletely reversed, that is, longer the molecule is, higher 
the conductance is. This special voltage is 2 V, where 
the order of current values is PA(5)DTs, PA(4)DTs, 
PA(3)DTs, and PA(2)DTs, in a descending order. 

Transmission spectrum of the biased molecule is an 
effective tool to understand the I-V characteristics. It 
denotes the position of HOMO and LUMO peaks and the 
related distance between the Fermi energy and HOMO 
or LUMO, so one can determine the point where the first 
steep increase of current comes out. 

As we know that the small HOMO-LUMO gap and 
the short distance between the LUMO (or HOMO) and 
the Fermi energy of the electrodes favor the higher con- 
ductivity. In Fig.2 HOMO-LUMO gaps of (n)TlDTs 
are very small and decrease obviously while number n 



increases. The peaks in the transmission spectrum or 
DOS spectrum illustrate the expanded molecular levels 
when the molecule is coupled with metal leads. Every 
expanded peak has tails overlapped each other which 
form the near-resonance channel. If HOMO peak and 
LUMO peak is very close, their strongly overlapped tail 
of transmission causes a big conduction value. Compared 
with other molecules, (n)TlDTs have very small HOMO- 
LUMO gaps with a great overlap in the low bias regime. 
This is the reason why (n)TlDTs have larger currents 
than other molecules under low bias. Furthermore, when 
number n increases, the HOMO-LUMO gap of (n)TlDTs 
decreases obviously. Hence, the overlap becomes stronger 
and the transmission in the low bias regime rises to in- 
duce an enhance of current. We can see clearly that the 
non-classical length dependence is just conduced by the 
drop of HOMO-LUMO gap and the distance between the 
HOMO and Fermi energy with increase of n number. 

Unlike (n)TlDTs, the HOMO-LUMO gap of 
PPh(n)DTs increases with the number of benzene 
rings, which indicates an ascending of the conductance 
gap with n number. LUMO, which is nearer Fermi 
energy than HOMO, dominates the position where 
the first resonance tunnelling happens and where the 
steep increase of current appears. The HOMO position 
of PPh(n)DTs is independent of n number. On the 
contrary, the HOMO-LUMO gap of PA(n)DTs decrease 
and the positions of HOMO and LUMO both approach 
to Fermi level when n number enhanced, which is similar 
to (n)TlDTs. Hence, the steep increase of current of 
long PA(n)DTs with small HOMO-LUMO gap will come 
out earlier than the short ones because its HOMO or 
LUMO resonance enters the bias window more rapidly 
than others. Under a little higher bias, a long PA(n)DT, 
e.g. PA(5)DT, throws out the resonance tunnelling and 
carries high current, meanwhile, the short molecule, e.g. 
PA(2)DT is still in the tunnel regime with low trans- 
mission coefficient. Therefor, the current of PA(5)DT is 
higher than that of PA(2)DT in the high bias regime. 
Because of such descending of HOMO-LUMO gap with n 
number, the I-V curves of PP(n)DTs will intercross, it is 
the reason why an irregular conductance oscillation and 
even the classical-reversed length dependence emerges in 
the high bias regime. 

At the same time, we notice that, in low bias regime, 
both PPh(n)DTs and PA(n)DTs have small current, 
which decreases exponentially with the rise of n number, 
because the low bias range is at the non-resonance tun- 
nelling regime between HOMO and LUMO. In Fig.2(b) 
the contribution from the tail of LUMO to conduction 
becomes smaller while n number increases. In Fig. 2(c) 
the situation is in reverse. The HOMO-LUMO gaps 
of PPh(n)DTs and PA(n)DTs are bigger than that of 
(n)TlDTs, which determines that the conductances of 
PPh(n)DTs and PA(n)DTs in low bias regime are mainly 
determined by the molecule length and thus exhibit the 
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facilities. 
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FIG. 3: Molecule bridge with the rectifying property. The 
resistances of PA(2)DT and PA(5)DT are controlled by the 
input voltage Vi n , which induce the sign variation in V a t- 



length dependence consistent with Magoga's exponential 
law. 

The intercross of I-V curves of PA(5)DT and PA(2)DT 
can be employed for designing the molecular bridge as 
shown in Fig. 3. V g denotes the intercross point. Be- 
low the point, the resistance of PA(5)DT is larger than 
PA(2)DT's, while above V g the situation will be in re- 
verse. When the input voltage Vi n exceeds the intercross 
point, V a b changes from the positive to the negative one. 

In a summary, we employ the first-principles method 
based on DFT combined with the nonequilibrium Green's 
Function theory to study the transport property of 
(n)TlDTs, PPh(n)DTs and PA(n)DTs. Conductance 
of the three oligomers shows the different length depen- 
dence of conduction. PPh(n)DTs are usualy like as clas- 
sical wires, but (n)TlDTs and PA(n)DTs demonstrate a 
crossover between the quantum conductance length de- 
pendence and the classical one controlled by the applied 
bias. The root reason is due to the fact that with the 
increase in the ring number in the oligomer compounds, 
HOMO and LUMO of the three oligomers approach to 
the different directions with respect to Fermi energy. The 
possible application of the length dependence of molecule 
conductance is proposed at last. 
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